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We propose an exactly solvable model for the two state curve crossing problems. Our model assumes 
the coupling to be a delta function. It is used to calculate the effect of curve crossing on electronic 
absorption spectrum and resonance Raman excitation profile. 
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INTRODUCTION 

Nonadiabatic transition due to potential curve (or sur- 
face) crossing is one of the most important mechanisms 
to effectively induce electronic transitions in collisions . 
This is a very interdisciplinary concept and appears in 
various fields of physics and chemistry, even in biology 
0. The most typical examples are, of course, a vari- 
ety of atomic and molecular processes such as atomic 
and molecular collisions, chemical reactions and molecu- 
lar spectroscopic processes. Attempts have been made to 
classify organic reactions by curve crossing diagrams Q • 
Other examples are nuclear collisions and reactions in nu- 
clear physics jj,|5j, dynamic processes on solid surfaces, 
energy relaxation and phase transitions in condensed 
phase physics 0, 0Land electron and proton transfer 
processes in biology [a, |fj . Two state curve crossing can 
be classified into the following two cases according to the 
crossing scheme: (1) Landau-Zener (L.Z) case, in which 
the two diabatic potential curves have the same signs for 
the slopes and (2) non-adiabatic tunnelling (N.T) case, 
in which the diabatic curves have the opposite sign for 
slopes. There is also a non-crossing non-adiabatic transi- 
tion called the Rosen-Zener-Demkov type 0, 0] , m which 
two adiabatic potentials are in near resonance at large 
R. The theory of non-adiabatic transitions dates back to 
1932, when the pioneering works for curve-crossing and 
non-crossing were published by Landau [1(1 ] , Zener [llj 
and Stueckelberg 0] and by Rosen and Zener re- 
spectively. Since then numerous papers by many authors 
have been devoted to these subjects, especially to curve 
crossing problems. Two categories might be classified for 
finding an exact analytical solution of the curve crossing 
problem. The first is that an exact analytical solution 
can be obtained for the whole region of the variable (say 
x here, see in the next section). For example, Osherov 
and Voronin solved the case where two diabatic poten- 
tials are constant with exponential coupling |14| . C. Zhu 
solved the case where two diabat ic p otentials are expo- 
nential with exponential coupling [l5j . The second is that 
an exact analytical solution is not possible for the whole 
region of the variable, but is possible for the asymptotic 
region. Then, physical quantities such as eigenvalues, 
scattering matrices can still be solved in an exact analyt- 
ical form, providing that the connection problem of the 
asymptotic solution is known. The Stokes phenomenon 



flif of asymptotic solution of the ordinary differential 
equation provides a powerful tool to deal with these kinds 
of problems [ij, Il8l Il9j | . Generalizing the real variable to 
the complex variable and tracing the asymptotic solution 
around the complex plane, the connection matrix which 
connects the asymptotic solution in the complex plane 
can be expressed in terms of Stokes constants. Recent 
work by Zhu and Nakamura found an exact analyti- 
cal solution of the Stokes constants for the second-order 
ordinary differential equation with the coefficient func- 
tion as the fourth-order polynomial. In this way, exact 
analytical solutions of scattering matrices were obtained 
for the two state linear curve crossing problem with con- 
stant coupling [2c| . 



THE MODEL 

We consider two diabatic curves, crossing each other. 
There is a coupling between the two curves, which causes 
transitions from one curve to another. This transition 
would occur in the vicinity of the crossing point. In par- 
ticular, it will occur in a narrow range of x, given by 



\V 1 (x)~V 2 (x)\ ~ \V(x c )\ 



(1) 



where x denotes the nuclear coordinate and x c is the 
crossing point. V\ and V2 are the diabatic curves and 
V represent the coupling between them. Therefore it is 
interesting to analyze a model, where coupling is localized 
in space near x c . Thus we put 



V(x)=K S(x-x c ), 



(2) 



here Kq is a constant. This model has the advantage 
that it can be exactly solved (delta function sink model 
for electronic relaxation in solution is known in litera- 
ture [2l|h One can roughly estimate the value of Kq 
as follows. If the two curves are approximated as lin- 
ear near the vicinity of crossing, then the energy sep- 
aration between the two is w [j§r(^i — Vi)] (x — x c ). 
This will be of the order of AV(x c ), if (x — x c ) takes 
the value ~ 7 „ A , VI ' x , c r \\ — . So it is reasonable to put 



Kn = 



2V(x c ) 



(&(Vi-v a )) 



x V{x c ). 
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EXACT ANALYTICAL SOLUTION 

We start with a particle moving on any of the two dia- 
batic curves. The problem is to calculate the probability 
that the particle will still be on any one of the diabatic 
curves after a time t. We write the probability amplitude 



(3) 



where ij}\{x, t) and ip2(x, t) are the probability amplitude 
for the two states. ^(x, t) obey the time dependent Schro 
dinger equation (we take h = 1 here and in subsequent 
calculations) 



dV(x,t) „ T/ , 
i—±±-L = HHf(x,t). 



H is defined by 



H 



Hxix) V(x) 
V(x) H 2 (x) J ' 



where Hi(x) is 



1 d 2 

H l (x) = - — — + V l {x). 
2m ox 



(4) 



(5) 



(6) 



We find it convenient to define the half Fourier Transform 
of by 



(7) 



Jo 

Half Fourier transformation of Eq. Q leads to 



Mo) 
MO) 



(8) 



This may be written as 

= iG(u;)#(0). 



(9) 



G(w) is defined by (ui — H) G{uS) = I. In the position 
representation, the above equation may be written as 



V(x,uj)=i G{x,x Q ]u))^{xQ,uj)dx , (10) 

J — oo 

where G{x,xq;uj) is 

G{x,x -lo) = (xKuj-Hy^Xo). (11) 

Writing 

G(x,xo; U )=(^ x '^ u \ ^f' x ^\) (12) 
\G 21 {x,x ;w) G 22 (x,x n -uj) ) 

and using the partitioning technique [22j we can write 

G n (x,x ;cj) = (x\[u-H 1 -V(u-H 2 )- 1 V]- 1 \xo). (13) 
The above equation is true for any general V. This 
expression simplify considerably if V is a delta function 
located at x c . In that case V may be written as V = 
K S = K Q \x c ) (x c \. Then 

G u (x,xo;m) = (x\[u - H i - KIG 2 ] {x c ,x c ]uj)SY~ 1 \xv) , 

(14) 

where 

G§(a;,io;w) = <x|(^ - ^J-^xq), (15) 

and corresponds to propagation of the particle starting 
at xq on the second diabatic curve, in the absence of 
coupling to the first diabatic curve. Now we use the 
operator identity 



(u-Hi- KlGl{x c ,x c ;u)S)- x = (w - Hi)" 1 + (w - fliJ-^GjftsBciEc; - Hi — K$G%(x c ,x c ;u)S}-\ (16) 



Inserting the resolution of identity 7 = \y) (y\ in gives 

the second term of the above equation, we arrive at 

G n (x,xo;u>) =G 1 (x,x - 1 uj) + K$G%(x,x c ;u) (17) W),- ncA 

*(%(z ,x eiU )Gu{xa,xo',u). (18) x[l-KtG%x c ,x c -,u)Gl{x M u)}-\ 

(20) 

Putting ,x = x c in Eq. (|17(l and solving for Gn(x c ,xo;u>) This when substituted back into Eq. I|17|) gives 
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KjGl{x,x c :uj)G%(x c ,x c ;uj)G<i(x c ,x ;uj) 
1 - KlG\{x c ,x c ;u)Gl(x c ,x c ;u) 



(21) 



r 



Using the same procedure one can get 



G 12 (x,x ;uj) 



1~K^G1(x c ,x c -uj)G^(x c ,x c ;uj) ' 



(22) 



Similarly one can derive expressions for G 22 (x, xq;uj) and 
G 21 (x,xq;uj). Using these expressions for the Green's 
function in Eq. © we can calculate ^(ui) explicitly. 
The expressions that we have obtained for ^(lo) are quite 
general and are valid for any Vi(x) and V^x). How- 
ever, their utility is limited by the fact that one must 
know G\(x, xq; u>) and G 2 (x,xq;u>). It is possible to find 
G®(x, Xq \ uj) only in a few limited cases and the harmonic 
oscillator is one of them. 



ELECTRONIC ABSORPTION SPECTRA AND 
RESONANCE RAMAN EXCITATION PROFILE 

In this section we apply the method to the problem 
involving harmonic potentials. We consider a system of 
three potential energy curves, ground electronic state and 
two 'crossing' excited electronic states (electronic transi- 
tion to one of them is assumed to be dipole forbidden and 
while it is allowed to the other) |23j, |24| . We calculate 
the effect of 'crossing' on electronic absorption spectra 
and on resonance Raman excitation profile. The propa- 



J 



gating wave functions on the excited state potential en- 
ergy curves are given by solution of the time dependent 
Schrodinger equation 




Allowed 



Configuraiional Coordinate 

FIG. 1: Schematic diabatic potential energy curves that illus- 
trate the model. The forbidden state is labeled F; the allowed 
state is labeled A. 



d ( vr b (M) 



dt V M ib (x,t) 



H vibiel (x) Vi 2 (x) 
V 21 (x) H vibte2 (x) 



r 



(23) 



In the above equation H vi t y ei(x) and H vi i, ye2 (x) describes 
the vibrational motion of the system in the first elec- 
tronic excited state (allowed) and second electronic ex- 
cited state (forbidden) respectively 



vib,el 



(x) 



l d 2 

2m dx 2 



+ \mu) 2 A (x - a) 2 (24) 



and 



H, 



vib,e2 



(x)=- 



J_d*_ 

2m dx 2 



+ -muj 2 F {x - b) 2 . (25) 



In the above m is the oscillator's mass, u>a and lof are 
the vibrational frequencies on the allowed and forbidden 
states and x is the vibrational coordinate. Shifts of the vi- 
brational coordinate minimum upon excitation are given 



by a and 6, and V\ 2 (V 2 i) represent coupling between the 
two harmonic potentials which is taken to be 



V 21 (x)=V 12 (x)=K S{x-x c ), 



(26) 



where K$ represent the strength of the coupling. The 
intensity of electronic absorption spectra is given by |24l 

E3 

I A (u) oc Relf^ dx n dx ^T b * (*) 

iG{x,x Q ;w + iT)*? h {x Q % (27) 

where 

G(x,x ;uJ + iT) = (x\[(LL> Q /2+uj-uj eg )+iT - e ] _1 |xo} . 

(28) 
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FIG. 2: Calculated electronic absorption spectra with cou- 
pling (solid line) and without coupling (dashed line). Here 
the values for the simulations are ujo = to a = luf ~ 400 cm -1 , 
T = 100 cm" 1 , s A = 10700 cm" 1 , e f = 11000 cm" 1 , 
K = 5.9643 x 10" 15 erg.A, and x c = 0.0895 A. 

and 

_ / H vi b,el(x) K \x c )(x c \ \ , . 

" mb > e - [k \x c )(x c \ H vibte2 (x) J [ZJ) 

Here, T is a phenomenological damping constant which 
account for the life time effects. ^ lb (x,0) is given by 

*f(x,Q)=( X f ) \ (30) 

where \i ( x ) is the ground vibrational state of the ground 
electronic state, loq is the vibrational frequency on the 
ground electronic state, ea is the energy difference be- 
tween the excited (allowed) and ground electronic state, 
and for the forbidden electronic state it's value is £p. 
Similarly resonance Raman scattering intensity can be 
expressed in terms of Green's function and is given by 

I R (u)<x \^ oo dx^ oo dx^f*{x,G) 

iG(x,x ;Lo + iT)^ b (x ,0)\ 2 . (31) 

In the above (m, 0) is given by 

*f(x,0)= ( Xf jf ] ) , (32) 

where Xf( x ) is the final vibrational state of the ground 
electronic state. As G®(x, xq;ui) for the harmonic po- 
tential is known j2^|, we can calculate G(x, Xo', Wc 
use Eq. i|31[l to calculate the effect of curve crossing on 
resonance Raman excitation profile. 



FIG. 3: Calculated resonance Raman excitation profile for ex- 
citation from the ground vibrational state to the first excited 
vibrational state, with coupling (solid line) and without cou- 
pling (dashed line). Here the values for the simulations are 
loo = to A = lof = 400 cm,- 1 , F = 100 cm -1 , e A = 10700 cm -1 , 
ef = 11000 cm -1 , K = 5.9643 x 10 -15 erg.A, and x c = 
0.0895 A. 



Results using the model 

In the following we give results for the effect of curve 
crossing on electronic absorption spectrum and reso- 
nance Raman excitation profile in the case where one 
dipole allowed electronic state crosses with a dipole for- 
bidden electronic state as in Fig. 2] As in |2J], the 
ground state curve is taken to be a harmonic poten- 
tial energy curve with its minimum at zero. The curve 
is constructed to be representative of the potential en- 
ergy along a metal-ligand stretching coordinate. We take 
the mass as 35.4 amu and the vibrational wavenumbcr 
as 400 cm -1 [24| for the ground state. The first dia- 
batic excited state potential energy curve is displaced by 
0.197 A and is taken to have a vibrational wavenumbcr 
of 400 cm" 1 . Transition to this state is allowed. The 
minimum of the potential energy curve is taken to be 
above 10700 cm -1 of that of the ground state curve. The 
second diabatic excited state potential energy curve is 
taken to be an un-displaced excited state. On that po- 
tential energy curve, the vibration is taken to have same 
wavenumber of 400 cto -1 . Its minimum is 11000 cm" 1 
above that of the ground state curve. Transition to this 
state is assumed to be dipole forbidden. The two di- 
abatic curves cross at an energy of 11673.6 cm -1 with 
x c = 0.0895 A. Value of K we use in our calculation 
is K = 5.9643 x 10" 15 erg.A. The lifetime of both the 
excited states are taken to be 100 cm -1 . The calculated 
electronic absorption spectra is shown in Fig. [3 The 
profile shown by the dashed line is in the absence of any 
coupling to the second potential energy curve. The full 
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line has the effect of coupling in it. The calculated reso- 
nance Raman excitation profile is shown in Fig. [3] The 
profile shown by the full line is calculated for the cou- 
pled potential energy curves. The profile shown by the 
dashed line is calculated for the uncoupled potential en- 
ergy curves. It is seen that curve crossing effect can alter 
the absorption and Raman excitation profile significantly. 
However it is the Raman excitation profile that is more 
effected. Even though our calculations are for a simple 
model (which has the convenience that everything can 
be found analytically), similar conclusions are valid for 
more complicated couplings [24| . 

CONCLUSIONS 

We have proposed a simple model for the two state curve 
crossing problem which can be analytically solved. Our 
solution is quite general and is valid for any potentials 
for which Green's functions for the motion in the absence 
of coupling is known. The model is used to calculate the 
effect of curve crossing on electronic absorption spectrum 
and on resonance Raman excitation profile. We find that 
Raman excitation profile is affected much more by the 
crossing, than the electronic absorption spectrum. 
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